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Abstract 

The low temperature configurational data obtained 
by Betts et al (1974) for the spin 1/2 Ising model on the 
hydrogen peroxide lattice are used to derive series in 
z= exp ( - 2J/kT), along the coexistence curve for the 
magnetization, M, and its first five derivatives, a M/ ou + 
where pw = exp ( - 2mH/kT). Series in wu along the critical 
isotherm are also derived for M and its first five deriva- 
tives, a M/az. Critical exponents and critical amplitudes 
are estimated by Padé approximant analysis. fecece £Otera 
few estimates which are not in good agreement with scaling, 
the majority of the estimated critical exponents support 


scaling theory. 
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CHAPTER 1 


Introduction 


C FaPTARD 
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Phase transitions are common enough occurrences 
that the layman is well aware of them. Yet his curiosity 
does not extend beyond a mere acknowledgement of the fact. 
The physicist however is impelled by a strong curiosity to 
get to the bottom of things. We know gases can undergo 
condensation to liquids under the right conditions; solids 
can be melted; and the components of a mixture in solution 
can be made to separate. Transitions of such assemblies 
belong to a special category which may be appropriately 
labelled as cooperative phenomena, or critical phenomena. 

We shall use both terms interchangeably. We use the des- 
cription "cooperative" because a detailed investigation 
reveals that the peculiar behaviour of such assemblies result 
from the interactions between individual microscopic systems. 
These interactions cause the cooperation of large numbers of 
systems to give rise to the phenomena listed above. In many 
of such phenomena, singularities or discontinuities are 
encountered in the thermodynamic description of these assem- 
blies, and so an alternative term for such phenomena has 

been "critical phenomena". 

The classical era of critical phenomena began with 
Andrews' 1869 publication of experiments on the critical 
point of carbon dioxide. A classic theoretical attempt to 
describe the critical region was made by van der Waals three 
years later, in his Ph.D. dissertation, "On the continuity 


of the liquid and gaseous states." To this day this des- 
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Cription is fairly accurate for temperatures that are not 
too close to Toe the critical temperature. Unfortunately 
there is a certain shortcoming in the van der Waals' 
theory. For all subcritical isotherms (T < a) there 
exists a region where the slope 93V/0P is positive, corres- 
ponding to the unphysical occurrence of a negative iso- 


thermal compressibility, K,, = -v"* (8V/3P) 9. The theory 


al 
was saved by the ad hoc construction named after its famous 
proposer, James Clerk Maxwell. For further details of the 
Maxwell construction one should refer to Huang (1963). 
Following the spirit of van der Waals' work, the 
next development of importance to the understanding of 
critical phenomena came from theoretical efforts to under- 
stand magnetic transitions. In the early twentieth century, 
important pioneering work on experimental magnetism was done 
by Curie, Hopkinson and others. Of interest to us are the 
physical properties of ferromagnets. At temperatures lower 
than a certain critical temperature known as the Curie tem- 
PeTacure, Tae a ferromagnet possesses a spontaneous magneti- 
zation. For temperatures above Toe the spontaneous magneti- 
zation vanishes, but the paramagnetic susceptibility becomes 
THLInL ce at Tor while the specific heat exhibits a \-type 
singularity. To explain and consolidate the experimental 
findings, a number of important theories were advanced. 
Pierre Weiss (1907) proposed a phenomenological theory of 
ferromagnetism in which he suggested the spins interact with 


one another through an effective molecular field. Moreover, 
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this field is proportional to the average magnetization. 

Later modifications of this idea of interacting 
Magnetic moments assume the constituent magnetic moments 
to be localized on fixed lattice sites, with pairwise 
interaction Of the “spins.*, In-parcticular, *the Ising moder 
has been exceptionally successful in explaining magnetic 
cooperative phenomena. The model was suggested by Wilhelm 
henz to his student EB. [sing in the ‘early 1920's; which is 
why the model is occasionally referred to as the Lenz- 
Ising model. The spins are assumed capable of two orienta- 
tions: "up" or "down"; with parallel alignment of neigh- 
bouring spins energetically favoured. 

Ising (1925) succeeded in solving the model with 
nearest neighbour interactions only for the case of a linear 
chain (i.e. a one-dimensional lattice with coordination 
number 2). The solution did not display a phase transition 
at finite temperatures greater than zero. Peierls (1936) 
demonstrated that a two-dimensional or three-dimensional 
model would exhibit ferromagnetism. This was substantiated 
in part by a major breakthrough by Onsager (1944) when he 
obtained the exact solution of the two-dimensional Ising 
model in the absence of a magnetic field. 

At the time of writing there is as yet no exact 
solution to the three-dimensional Ising model. Fortunately, 
approximation methods are available in the form of series 


expansions of the partition function at low and high tem- 
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peratures (Domb, 1960). More recent reviews can be found 

in Fisher (1967), Stephenson (1971), and Domb and Green 
(1974). With the aid of high temperature series expansions, 
the critical temperatures and the critical index, y, of 
PNeginlelLale susceptibility for thestsccc., bscvc., S.c. and 
diamond lattices (coordination numbers 12, 8, 6 and 4 res- 
pectively) have been investigated and are well known for 
some time now (Domb and Sykes, 1957; Essam and Sykes, 1963). 
However investigations on ae low temperature side have been 
quite scarce until recently (Essam and Hunter, 1968; Betts 
and Filipow, 1973; Betts and Chan, 1974). 

The general theory of low temperature series expan- 
sions have been laid out in a series of papers: Sykes, 
Essam and Gaunt, 1965 (paper I); Sykes, Essam, Gaunt and 
Hunter, 1973 (paper 11); and Sykes, Gaunt, &ssam, Mattingly 
and Elliott, 1973 (paper III). Now the hydrogen peroxide 
lattice (Heesch and Laves, 1933; and Wells, 1954), is the 
three dimensional lattice with the smallest coordination 
number (q = 3). Using it as a lattice for extending the 
three-dimensional model to the lowest limit of coordination 
number possible for a three dimensional lattice, Leu, Betts 
and Elliott (1969) investigated the critical properties on 
the high temperature side (T > TQ): The high field (low 
temperature) series expansions for the Ising model on the 
hydrogen peroxide lattice were derived in a later paper by 


Betts, Elliott and Sykes (1974). We shall make use of the 
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configurational data obtained by Betts et al (1974) to 
investigate the following. For T < Toe we shall study the 
critical behaviour of the magnetization M, and its first 
five derivatives with respect to the field variable, 

u = exp (-2mH/kT) on the coexistence curve (H = 0). We 
shall also derive series expansions in uw along the critical 
isotherm (T = Toe for M and its first five derivatives 
with respect to the temperature variable, z = exp(-2J/kT). 
The techniques of Padé approximant analysis will be applied 
to the series so obtained to estimate critical exponents 


and critical amplitudes. 


CHAPTER 2 


Theomy.of the critical region. 
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2.1 Critical exponents 

In the study of critical phenomena it has been 
found that most of the thermodynamic functions of interest, 
such as the specific heat, the magnetization and its higher 
derivatives, usually vary with temperature or some other 


independent thermodynamic variable according to the form: 
f(x) Af (x-x,) /x}” +less singular terms, X>X 


A'{ (x -x) /x}"+less Singular terms, X<K Oy (Fi) 
where A, or A', a numerical constant, is called the critical 
amplitude, x is a general, independent thermodynamic variable, 
x is the value of the independent thermodynamic variable at 
the critical point, and f(x) is a general thermodynamic 
function of interest. Then the numerical constant, i’, or ir', 
is defined as the critical exponent characterizing the asym- 
BrotLe behaviour of the thermodynamic function, E(x), at the 
critical point. In view of the less singular terms in (2.1), 
the critical exponents, or \', are sometimes defined more 


rigorously as: 
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We look next at a few selected critical exponents 


and amplitudes. We define t= (T-T )/T. , where T is the 
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temperature, and Tt. the temperature at the critical point. 
Since we are interested in reyions close to the critical 
point we ignore less singular terms in the definitions 

that follow. For a magnetic system the specific heat at 
constant magnetic field, H = 0, has the following asymptotic 


functional form near T = T 
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The dimensionless quantities A and A‘ are the critical amp- 
litudes. The critical exponents a and a' characterize the 
behaviour of the specific heat at constant field when the 
critical point is approached from above and below To res- 
pectively. For the three dimensional Ising model it is 


believed that (Sykes, Hunter, McKenzie and Heap, 1972) 
Cee L/S (251) 


The temperature dependence of the reduced zero 
field magnetization, M(t) near the critical point is 


described by 
My (t) = My (t)/ ty (-1) ~ B (-t)?, <0 a-)) 


(Jn is used here to denote the magnetization in the ordinary 


sense.) 
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We note that B, the critical amplitude, is again dimension- 
less and 8 is the critical exponent characterizing the 
behaviour of the zero field magnetization in the neighbour- 
hood of the critical point. The estimates of 8 (Fisher, 


29Gs)) Jae atettesrange ).303°* Ge (0.318, The conjecttirs 1s 
fee (5 Po = 907312500 (2:6) 


The exponent 6 describes the variation of the 
reduced magnetization, M, with the reduced field (h = mH/KT .) 


om the Critical rsotherm t=0 (T=T : 
Im] ~ Blnj2/® (2.7) 


(k is Boltzmann's constant, and m is the magnetic moment 


per spin). Sometimes the above expression is written as 
h wD M (2.8) 


it is believed (Gaunt and Sykes, 1973) that for the three 


dimensional Ising model: 
O25 C229) 


Consider the H-M isotherms for a magnetic system. 


The slopes of these isotherms are proportional to the 
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to infinity as T approaches Toe So we define critical 
exponents, Y and y', and critical amplitudes, C and C' 


to describe this behaviour. 
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Xp is the isothermal susceptibility of a system of non- 


interacting magnetic moments: (paramagnet) at the critical 
point. It is believed that (Fisher, 1967; Sykes, Gaunt, 
Roberts and Wyles, 1972): 
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Zine Exponent Inequalities 
The work of Rushbrooke (1963), Essam and Fisher 


(1963), Fisher (1964, 1967, 1969), Griffiths (1965a,b, 1968) 
and Griffiths, Hurst and Sherman (1970), among others, 
have focused attention on a number of rigourous relations 
amor.g critical exponents expressed in the form of inequal- 
ities. 
The Rushbrooke Inequality 

For a magnetic system we may define specific heats 


at constant magnetization, C or at constant field, C 


M’ H* 
them: 
Oe ee aR 
C2 tT 35) Se ae (aac) (2.13) 
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M 
and: 
aS JE ae 
H H oT 
H 
The isothermal susceptibility is given by: 
2 
_ ( oM o 0G . 
ue a) ae (—>) (2.15) 
dk oH T 


while the adiabatic susceptibility is: 
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In the four expressions above we have: 
U = internal energy 
» = entropy: 
A = Helmholtz potential 
G = Gibbs potential 
BE? =enthalpy 


From thermodynamic considerations it can be shown that 


_ aM 
Che Cher Fea (31) (oe) 


Griffiths (1964) showed that (i) the Gibbs potential, G(T,H), 
is a concave function of both magnetic field and temperature, 
and (ii) the Helmholtz potential, A(T,M), is a concave 
function of the temperature and a convex function of the 
magnetization, for a magnetic system whose Hamiltonian is 


of the form: 
A = AL, = Het (2,18) 


cM is the magnetization operator. 

(The Ising model has a Hamiltonian of the form (2.18).) 
Using the convexity property of A(T,M) it can be shown that 
Cu must be positive. Then taking equation (2.17) it is 


found that 
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From (2.19) and the definitions of the critical exponents 


discussed in section 2.1, Rushbrooke (1963) showed that 


ging BePtrey' tLaea, (2.20) 


(2.20) is called the Rushbrooke inequality. 


The Griffiths Inequality 


The other inequality of interest to us is the 


Griffiths inequality: 
GPe+OR (Lip ee per’ 2% (221) 


We omit the proof here, and instead refer the reader to 


Griffiths (1969a, 1965b)'. 


rsnscospe fnotise aie pr ee ong bas non 2) me +1 
‘nisi 


tans beworde | (eeet) oxinondelebl esha nokthse at boeauoe 
a 


dj 


— 


(OS.8) | Se ye ast 


yditsupeni stoordiesA ett pefies at (0S,S) 
yi iisupent efjittis oft 
les - 
adj et ev of Saorstnl to vt feu: ant tenhito ait 


‘yd iLeupiens ort 19083 
(f8.8) 8 < Co 4-18 +O + 


od xyebser oft sete bseteni bae ,oyor Tooxgd ans = ow 


(deel ,e@aes) ads t?%, bx 


Le 


ORE Scaling Ideas 


An immediate question is, "When are the inequal- 
ities discussed in section 2.2 equalities?" Approaches 
investigating whether RASES inequalities are indeed 
equalities were first described by Widom (1965a,b), 
Patashinskii and Pokrovskii (1966), Domb and Hunter (1965) 
and Kadanofé (1966). Today these ideas are embodied in 
what has come to be Heese iced as scaling theory. 

It should be pointed out that standard scaling 
theory is essentially a theory for critical exponents. It 
has also come to be called the homogeneous function approach. 
At the heart of this approach is the assertion that the Gibbs 
potential, or free energy, is a generalized homogeneous 
function. Then from the form of the Gibbs potential and the 
properties of homogeneous functions it becomes apparent that 
all of the critical exponents can be expressed in terms of 
two parameters, if we consider a magnetic system where the 
Gibbs free energy is a function of the temperature and the 
applied field. Two critical exponents, 8 and 6, have 
already been introduced in the preceding section. Two 
methods will be used to show how the critical exponents of 
the higher derivatives of the magnetization are related to 
B and 6. 

We begin with the approach due to Widom(1965a,b). 
The basic postulate of static scaling theory tells us that 


the Gibbs free energy is a generalized homogeneous function 
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of the field, and the temperature. Without loss of 
generality, we shall work with the reduced temperature, 
LoS (T-TY/T Oe and the reduced field h = mH/KT .. Intro- 


ducing two general parameters, a, anda the Gibbs free 


te h’ 
energy may be written: 


GOVT ae ihe) = A eli h) C22.) 


The above is then differentiated with respect to h: 
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r 2 ateq t 


But the reduced magnetization is defined as: 


M(t,h) = - 9G(t,h)/dh (2e24) 


Hence (2.23) simplifies to: 


h M(A “t, A “h) = 2 M(t,h) (2.25) 


It now becomes apparent that there are 2 critical exponents 


associated with the behaviour of the magnetization near the 


critical point depending on whether 
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cea h > 0, equation (2.25) becomes: 


at a 


M(Omh) =A m(0,2."h) (2.26) 
-l/a, 
Secting AJ= h , one obtains: 
\i=a,)/a 
Meh) eh ee a) Oa 
eas 


Buc trom (257):; M(0,h)-= h 
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When *# “h/iA ~t = 0), equation (2.25) reduces to: 


M(t,0) = A M(r ee 10) (Aveo) 
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Now upon choosing i = (-1/t) , we get: 
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becomes such that: 


Mier Oy tee (2.30) 


(Note that less singular terms are neglected as before.) 


Thus B = (l-a,)/a, (2075.1) 


Solving for a,: 


a, = pot ter ene 
Seonving (2.27) Lor a: 
ay, = C7 Coe) (23.33) 
Substituting for ay and ay in equation (2.25) yields: 
eye te ee II Pe ey sd) 
With a choice of A = Us , (2.34) becomes: 
m(t,hy = §-¢m(s47%t, 5°) eesigl 


The above can be simplified into two convenient expressions 


by either choosing § = pees Go eraann ts © 
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Since we are interested in the low temperature (T < TQ) 


side of the critical point, equation (2.36) is modified to 


read 


M( th) = (-t)® ¢ (hv(-t)>) (2.37) 


1/6 


When a choice of § = h- in (2.35) is made: 


aU) seus Sesane = as (2.38) 
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The above is consistent with the form obtained by Betts and 


Filtpow ((1972:)i. 
A second and shorter approach is to say that in 


general, for a magnetic system, one may write: 
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Scaling theory tells us that instead of (2.39) any two of 


the variables h/m® ; m/t?, Hy tee may be used to describe 


the functional form of M or h. 


Thus one may write: 
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h/M = F, (m/e!) (2.40) 
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or 


h/m° 


tt 


£, (h/t”) : (2.41) 
or 


m/t? 


£ (h/t) : (2.42) 


Also scaling theory tells us that as h/t = .0, f, > constant 
in (2.42) and, modifying the equation for t < 0, we find 


along the coexistence curve: 
My(t) ~ (-t)? (2.43) 


The isothermal susceptibility, Xqe is given by 


Xp = (dM/dh) ,,. From (2.25), on differentiating with respect 


to h keeping T (and hence t) constant: 
eee ene) (2.44) 
But by definition of the critical exponent jy: 
Pane (2.45) 
Setting h/t’ + 0 in (2.44), 


ante Se (2. 46) 


Comparing (2.45) and (2.46): 


re re 
batt ow 0. 4 sok seddetps ¢ higeieee 1 Die (ap. 8 
ae ‘rita 


(Et ss} aa) 7 Gin 


yd movip Bh ok ich. Haein Nea RS PT 7 
sosgeskt dtiw ondtsetsterottib: no , (AoyS)" mnt vig (HENNEY = ok i . 
OUP 


:Yasdanos (4 sonet bos) Topeiqesd a : 
— 


(pe oh PMN cl tn 


Y = 4-68 (2.47) 


In three-dimensions no critical exponents have been estab- 
lished exactly. The work of Sykes, Gaunt, Roberts and 
Wyles (1972) however, strongly indicates that the suscepti- 
bility exponent, y, is very likely 5/4. We shall assume 


y = 5/4. 


Thus: A-B = 5/4 (2.48) 


On differentiating (2.22) twice with respect to t, and 


using (2.14) and (2.3), it can be shown that: 


a+ 28 + y = 2 (2.49) 


As mentioned earlier in section 2.1, there is every indica- 
tion from Sykes, Hunter, McKenzie and Heap (1972), that a 
is probably 1/8. We shall assume a = 1/8. Then using (2.4 


and (2.48), we have: 


Beas S/P65e and Av =—25/16 (25550!) 


The next step is to compute the critical exponents for the 


higher derivatives of the magnetization. Two cases are 


considered. In one case the critical point is approached 
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Figure 2.1 


Paths along which the critical point is 


approached. 
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along the critical isotherm, In the second case it is 


approached along the coexistence curve. See figure 2.1. 


Case (1) critical isotherm, 
To examine the temperature derivatives of the mag- 
netization on the critical isotherm, it is more convenient 


to write (2.38) as: 
m(t,h) = hP/* g(tynt/4) (2.51) 
Ignoring less singular terms: 


(SEM 9s) Sedan ke ee (2.52) 


i.e. /(0°M/dt'),_4 ¥ hk (2.53) 


where €, = (2-B)/A. Table 2.1 gives the scaling 


Preedietions fore Veil ie pod. Os 
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Case (ii) coexistence curve 


Recall equation (2.37): 
nen oman o | 


Therefore, taking derivatives with respect to h, and 


omitting less singular terms: 
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ates? B-2A 
(an ag ae) (2.54) 


i.e. (@°M/ah’),_. v (-t) Ye (2.55) 


where 9 = §A=68. Table 2.) gives the scaling predictions 
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CHAPTER 3 
The three-dimensional Ising model on the 


hydrogen peroxide lattice 
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cee Configurational Energy 


We present here a brief outline of a method due to 
Sykes, Essam and Gaunt (1965) for deriving the free energy 
of the three dimensional spin 1/2, Ising model. Let 2J 
be the energy gained if 2 first neighbour spins alter 
alignment from parallel to antiparallel orientation, H be 
the applied magnetic field, q be the coordination number 
(3, in the case of the hydrogen peroxide lattice), and m 
be the magnetic moment per spin. At absolute zero there 
is complete order and all the spins point one way, thus 
giving rise to a spontaneous magnetization. The energy of 


the ground state for N sites is 


-N(1/2 qJ + mH) (3.1) 


As the temperature rises, thermal disorder occurs. 


The overturning of s spins with r first neighbour bonds will 


lead to an increase in energy of 


2(qs - 2r)J + 2msH (Gia) 


We define: 


Z exp (-2J/kT) 


i} 


7 exp (-2mH/kT) (a0) 


where k is Boltzmann's constant, and T is the temperature 


in degrees Kelvin. 
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The probability for the thermal perturbation is 
given by the Boltzmann factor: 


zas-2r hs 


exp{[-2(qs-2r) J ~ 2smH] /kT} — (3.4) 


Let Ay be the partition function for N sites, and A he 
the partition function per site. For large N, Ay av ies 

At low temperatures A, which involves the sum of product 
of factors like those appearing in (3.4), can be expanded 
in terms of successive deviations (corresponding to more 
and more overturned spins) from the lowest energy state, 

as Domb (1960) has pointed. So the free energy per spin 


may be written as (log implies taking the natural logarithm): 


¥y = -5 qd - mH - kT log A(u,z) (335) 


Once the partition function is known, the formalism of 
statistical mechanics enables us to calculate the various 


thermodynamic functions of interest. 
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3.2 The Low Temperature, Hich Field Polynomials 


We are now faced with the problem of obtaining a 
series expansion for the natural logarithm of the Partition 
function per site. To derive this series, perturbations 
of the ordered state must be studied. In general, the 
number of perturbations of a given configuration on a 
large lattice of N sites will be some polynomial in N. 


Domb (1960) showed that Mfu, 2) can be expanded in the form: 


Mu,z) = 1 + uF, (N,2) + w°F, (N,2)+...4uSP_(N,z)... 
(3.6) 
Each F(N,2) is a polynomial which can be expanded in 
powers of N. Domb also showed that if instead an expansion 


of logA(u,z) is preferred, then: 
Log (Zz) = uL, (2) + UL, (2) toot WPL. (2) +... Pere) 


where L. (2) is but the coefficient of N in FS (N,z). L, (2) 
is itself a polynomial in z. The set of polynomials, 
{Lo (z)}, is the set of low temperature (high field) poly- 
nomials. 

To specify a given L. (2) we require the contribu- 
tions from perturbations with s overturned spins. The s 
spins may occur with zero first-neighbour bond, 1 first 
neighbour bond, 2 first neighbour bonds, ... and so on. 


Accordingly, following Sykes, Essam and Gaunt (1965), we 


write: 
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58 (s-1) 
ie ea 5) [s;x] 295°2 (3.8) 
=0 
The complete high field polynomials for the hydrogen 
peroxide lattice for s < 24 were obtained by Betts, 
Elliott and Sykes (1974), and are presented in appendix 
2 for reference. 

Take the hydrogen peroxide lattice. Consider a 
single Site of the lattice. Then if the three nearest 
neighbour sites are taken as vertices of a triangle it 
is possible, by omitting the site we started with, to 
derive a sublattice of second neighbour vertices. The 
sublattice so obtained is called the hypertriangular 
lattice by Betts, Elliott and Sykes (1974). The hyper- 
triangular lattice derived from the H,0, lattice is 
analogous to the triangular lattice derived from the 
honeycomb lattice. In appendix 1 the hypertriangular 
lattice graphs and low temperature lattice constants needed 
for obtaining Fo(s < 12) are given. 

On examining the hydrogen peroxide lattice one 
sees that it may be easy to determine visually the number 
of ways in which a polygon of ten edges may be embedded in 
the lattice. However as polygons of more edges are con- 
sidered, the task of deciding how many ways a polygon may 
be embedded in the hydrogen peroxide lattice becomes more 
difficult. To handle this difficulty sophisticated book- 
keeping has to be developed. To see how this is done, one 


should refer to the work of Leu, Betts and Elliott (1969), 
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and Betts, Elliott and Sykes (1974). 
Armed with the low temperature (high field) 
polynomials we are now ready to derive various thermo- 


dynamic functions of interest from the free energy given 


shgy (Se is 


Cys The Ma netization and its Derivatives 
By definition the magnetization per spin is given 


bys 
M = ~(39/9H) 9 
We shall work with the reduced free energy per spin: 


Go= ¥/kt 


= -(q logz )/4-hT_/t- J Lo (z)u® (9) 
S 


Thus we have for the reduced magnetization per spin 


(h = mH/KT |) : 


= 
I 


= (9G/9dh),, 


— (06/9) 5 (dun/9h),, (3.10) 


ieee } s Lo (z)u 


1.66. BG 2) 


Among the possible paths one can take in approaching the 
Critical point, two are particularly convenjsent., 0nesis 
a path along the coexistence curve and the other is along 
the critical isotherm. Figure 2.1 shows these two paths. 
° LJ LJ s 2 b 
The estimate of the critical point z, 1s necessary efore 
one proceeds with the analysis of the Ising model on the 
hydrogen peroxide lattice. Leu, Betts and Elliott (1969) 
have estimated sie tats tanh (J/kT), from which they obtained 


2 = 0.317401 + 0.000010 for the hydrogen peroxide lattice. 
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They did this by analysing the high temperature expansion 


of the susceptibility. 
(a) on the critical isotherm. 


By setting z= Zo in (3.10) one obtains: 


SUC co Ae lal cree (ese (3212) 
Ss 
In general the &th derivative of the magnetization with 


respect to z can be written as (2 >0).: 


[a°M(u,z)/az") 2 -23u>[a°L (2) a2") elon 
Cc S} (e 


ImMenerore upon summarizing (3.12) and (3.13) 
23 
[a'M(u,z)/a2"1 5, = } aoe a (3.14) 


The summation is truncated at s = 23 because at the moment 
we know L, (2) only up tos; = 23. The series) (3214) for 
l= 0 to are given in Appendix 3. The coeflicients are 
quoted to only six figures because they are dependent on 
Zo which is known to only six figures. 
(b) on the coexistence curve. 
When (3.11) is differentiated with respect to iu, 


keeping z constant we have for the %£th derivative: 


23 : 
gem (u,z)/aue = -2 ) s@(s-1)...(s-2+1) Li (z) yu” x (3.15) 
s=k 
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Along the coexistence curve we set up =u. =1 


(zero field, h = 0), and write: 
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The reason for truncating the series in z is that we know 
L. (2) only up to s = 23. Each L, (2) when written out in 
ascending powers of z, begins with a leading term w(s,p)z?, 
where w(s,p) is the first non-zero numerical coefficient of 


ae 


in Li(z). When s is tabulated against p we notice that 
the trend indicates that with Lo,(2), p18, 18. So Lt ws 
sensible to retain terms in z up to a only. Hence the 
series infz in’(3.17) are truncated at zal The series for 


tee ou kb, 2.4.5 im 1(3,.17). care: given an, Appendix 4% 


Table’3.1.°' A table of p vs. s where p is the smallest 


power of z in L(2). 
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From (2.52) and (2.54), scaling theory predicts 


that: on the critical isothern, 
u x 43 
(O° °M/ot ae, he ee i 
and on the coexistence curve, 
Q te 7 oe -y! 
(0° M/oh Vea (-t) ‘% ; 


Since our series expansion variables are u and z, then, 


provided less singular terms are omitted: 


Re hee 
(Q°M/3z") _o 


(a"M/at"), 4 ~ (az/at)* ; 


eco) 
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Hence by studying the asymptotic behaviour of (0 M/dz ae 
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CHAPTER 4 


Analysis for Critical Exponents and Amplitudes 
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4012 Pade Approximants 


We have seen that the magnetization and its 
derivatives with respect to z and u are expressible in 
the form of series expansions. Since the high field poly- 
nomials L. (2) are known only to s = 23, there are unavoidable 
truncations of terms of the series for (a°M/az") and 
(Q°M/Iu") - Our problem then is to extract the Serre 
behaviour RA a given finite series. We hope that in spite 
of its finiteness it already contains the characteristics of 
the critical behaviour. This is not a vain hope as the 
success of several workers in this field has shown. 


co 


Consider a function f(x) = ) wx. Suppose we 
n=0 
truncate it at n = L. The Padé approximant to f(x), [N,D], 
(Baker, 1961; Gammel, Marshall and Morgan, 1963; Baker, 1965) 


is the ratio of 2 polynomials of orders N and D respectively, 


yates ND) 2) = —— eee (4.1) 


such that 


Lowe} vA 
WH ee px n 
om i n=0 a n=0 q,* oe) 


The coefficient do is chosen to be 1 without loss of 
generality. If each Ne is known (ne=-.0,,.., 2a.) sscnen 


the other coefficients can be determined provided N+ D < lL. 


sachin, imate 8s ~ <0 GB teleld ‘| 
bis ssn ay ao\M ED 162 etree afd %o winted Io enokd ae 
fesisrr5 sae dostaKe of Bt nets matdorg 05 * a a toe ‘ he 
etig2 ni add Sqod ov -e5irse Perec’ asvip 5 mot? suoiveris do 


i> esiveiveineyento sia entetnds ybseris $i ecanetint? ati 2 > 
= 


eid #6 sqné aisvs Jon ei eft .idotveded Inotstso: sian 
.twofle est cielo aids oi avexstow lpusvee Yo eaeont ya 


yr "Ry ? _ = i noitonut 6 18B£eno9 # _ iF 


-_~?¢ 


, 1G, Mi sie ot Pee ey em S6n¢ sit .d = a ts 4% © 


(2501 ,xeded 298% /akpxoM Ban Wédeiem ,lommgn L304 , x9468) 
\Ylevisesqass @ Bas W eYshao to aleimonylog ¢ t¢ obtss ad al 


(Ef. by = i) fay) vest 


38 


This is done in the following manner. From (4.2) we have: 


L De i N 
(wo + W)X +..-t Wy x ) (1 + q,* +...+ Gpx ) = PotPy xt. - -+Pyx 


By equating coefficients of the constant term up to x we 


find that: 
WO a 
and 
An bale Gy Py 
Wy ee q2 P2 
e + ° e ° = C) (4505) 
Wy Weed Newt Ou ctN Py 


where, if D < N, we get Ine’ Ape2' °° A = 0. 


By equating coefficients of ues to hey we have: 
ASR WN Nad ie hart |: 
NED a] Sean WN=D 0) 
: a ae ‘ 3 z = 0 (4.4) 
WN+D Wy+p-1 “N+D-2 Wn ap 


where if N < D we set W_,, Wy --- to zero. Equations 


(4.3) and (4.4) may be solved on the computer. 


7 
. / 4+ OR ,oxedw AF 

<0 = oP ‘* * Acag? ‘page: 282, 24 ‘ | | a . 
:ovad sw om, od ru. ‘0 fi 7eio0r7 teen oe 7 7 


— 


a.) o = 
> 


| ot rot oe Ww 3 ase ew od eee 


She 


Due to the Padé approximant being a rational 
function it is very suitable for applications to thermo- 
dynamic functions which have a power law singularity. 
For example we may consider 


ata nia py 7¥ 
X(T) = d Ane ACE a1) (4.5) 


Taking logarithmic derivative of the above yields: 
dlogx(T)/dT ~ y/(T_-T) (4.6) 


We see by the expression (4.6) that the Padé approximant to 
the logarithmic derivative of x(T) will have a simple pole 
at T= T. and that the negative of the residue of this Padé 
approximant is the critical exponent, y. 

Suppose the negative of the residue (i.e. y) is 
plotted on the ordinate axis, and the pole (an approximation 
Sor TO) on the abscissa axis. Each Padé approximant to (4.6) 
provides us with a point. A smooth curve can be drawn 
usually through the collection of points due to the set of 
Padé approximants, {[N,D]}, used to approximate (4.6). The 
best estimate of y is that ordinate on the curve for which 


the abscissa is the true Ln (assuming T. is known). 
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AyS2 Logarithmic Derivatives 


We now wish to apply the method mentioned in the 


preceding section to the series (2 = 0, 1, 2, 3, 4, 5): 


23 
Q Q iz (2) n 
fa ean pee ees eau (4.7) 
c n=0 
and 
as eres are 2 5 (2) on 
C 10} ear = n Cir. (4.8) 
c n=0 


Case (a). 
The most singular parts of (4.7) are expected to 


diverge as 


~€ 
(a°M/az") ~~ EQ-u) * (4.9) 
Cc 


When the logarithmic derivatives to (4.7) are formed we 


would expect them to diverge as 
Q Q 
(d/du) {log (ad M/az Ven a E,/(1-u) (4.10) 
¢ 


Each Padé approximant to (4.10) should yield a pole close 


to have 1, and a corresponding residue that is an estimate 
of Eps 
Figures 4.1 to 4.5 are plots of the negative of 
Ch) on 


the residues of the Padé approximants to d{log(ta u )}/dy 


n 


versus their corresponding poles closest to the physical 
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Singularity, Ho Ly tor = 0; 1, 25.4, 55. The results 
for £ = 3 are so scattered that the meaningless plot has 
been omitted. Note that enn -1/6, where 6 is defined 
in section 2.1. In table 4.1 estimates of Ey with their 
confidence limits, as obtained from the graphs, are listed. 
It can be seen from figure 4.1 that the points are 
crowded into a region well below lea Ll. “This makes: tor 
difficult estimation of Ey: Similarly the concentration 
of points above ue = 1 in figure 4.2 renders a precise 
estimation of ey impossible. The E,'s so estimated 
(2 = 0, 1, 2, 4, 5) are used as guides to the choices of 
powers made in a more precise technique that of neighbouring 
powers. 
Case (b). 
The analysis of (a°M/3u") proceeds in a similar 
fashion to that of (3'M/a2") ; the most Singular parts 
Cc 


of 4.8 are supposed to diverge as 


ee =o 7 Q 
(d°M/dy usu. ~C, (1 z/2.) (A511) 


When the logarithmic derivatives to (4.11) are computed the 


following divergent behaviour would be expected: 


(d/dz) (10g (2"M/ 3H) yay ~ Yo" / (2-2) (agi2) 
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A given Padé approximant to (4.12) may have a number of 
poles, but the pole of interest is that one closest to 
Z= 2Z_. The negative of the residue computed for this 
pole is plotted against the pole. Figures 4.6 to 4.11 are 
such graphs of the negative of the residues versus the poles 
of the Padé approximants to A{log (Zb*) 2") } /az for Pah 
Pear Opel 2, ogee, oe The ‘estimates.of a) obtained from 
these graphs are displayed in table 4.2, together with their 
confidence limits. 

In figure 4.6 the poles cluster well below 
Base 0.317401, making it difficult to pinpoint Yo precisely. 
It is readily seen that the scatter in the points of the 
graphs used for estimating Y4 and Ys are quite bad (figures 
4.10, 4.11). One possible reason for the erratic behaviour 
of the Padé approximants used for these estimations is the 
presence of non-physical singularities. These non-physical 
singularities prove troublesome when their locations in the 
complex plane of the independent variable (in this instance 
the complex z plane) are nearer to the origin than Zs OF 
are clustered very close to the physical singularity at 
Z=2z_. In a subsequent section we will look at transfor- 
mations that may facilitate analysis of series expansions 


with such problems. 
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Figure 4.1 Estimates of WH, and EQ from poles and 


residues of Padé approximants to (d/dy) log M(u,z_). 
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Figure 4.2 Estimates of u_ and ey from poles and residues 


of Padé approximants to des log Z [Bayes 5 
Zz 


du U OZ 


Pigure 4.3 
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Estimates of We and E> from poles and residues 
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Estimates of We and zl from poles and 


residues of Padé approximants to 
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Figure 4.5 Estimates of Ho and Eg from poles and 


residues of Padé approximants to 
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Figure 4.6 


Estimates of Z and Yo from poles and 


residues of Padé approximants to 


de | 109 mc}. 


‘3166 ‘370 


‘3174 


48 


Figure 4.7 
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Estimates of ZO and Yq from poles and 


residues of Padé approximants to 
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Figure 4.8 
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Estimates of Zo and Y5 from poles and 


residues of Padé approximants to 
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Figure 4.9 


Estimates of Z and y, from poles and 


residues of Padé approximants to 
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Figure 4.10 


Estimates of Zo and v4 from poles and 


residues of Padé approximants to 
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Figure 4.11 


Estimates of Z and Ye from poles and 


residues of Padé approximants to 
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were made using the technique of computing Padé approximants 


to the logarithmic derivatives of (4.8) and (4.7). We 


recall the divergent behaviour of the most singular parts 


ofe(4.8), and «4<67)= 


Oo Muy) 2) 
Oz 
Z=Z 
(e 
Q 
d M(u,Z) 
0 
Pp) 
4: w=y 


On raising (4.9) to the 


1/Yq1 we get 


3M (u,z) 


s2 


C 


and 


aM (u,z) 


Q 


1/Y9 
du | 


From tables 4.1 


eand Yo by the method 


Ne 


Call these estimates eo 


neighbouring powers, we 


=F 
B, (1-w) ‘ (4.9) 


Ci iy aeei: 
Clea a (4.11) 


power l1/e and (4.11) to the power 


fey 
1/e, 
os ey) /(1-u) (4.13) 
aye /On5 
shes) GRY BAS) (4.14) 


and 4.2 we have best estimates of 
of taking logarithmic derivatives. 
and Po: In the technique of taking 
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Q toaeer 
lo M(i, 2) /a2 lite for a series of values of Ey in the 
e) 
neighbourhood of Eo. For each Eo and each Padé approxi- 
mant to (4.13) we get a pole which is an estimate of Wor 


NES 
and a residue which is an estimate of ve 2 Lt, 2s 
then possible to plot a graph of the estimated Wu (IN,D]) 
versus ©) for each Padé approximant, [N,D]. The best 


estimated €, in this technique is that Ep which corresponds 


iv 
ie 8) we de: 
Similarly, for a fixed 2, we form Padé approximants, 
Q pa? tp 
t{N Yep leb to tho M(u,2Z)/on-] for a number of values of 


Yq in the neighbourhood of Vo. This time, for a chosen Von 
the simple pole of the Padé approximant to (4.14) yields an 


estimate of Zoe while the residue calculated at that pole is 
L} 


en 
an estimation of (CT) be It is again possible to plot for 


each Padé approximant, [N',D'] a graph of z({N',D']) versus 
Yo: The best value of te is that value which best reproduces 
ZN 0.317401. It is clear from an inspection of figures 
4.12 to 4.22, how we arrive at confidence limits in the 


estimation Of ©, “and Yo by this technique of neighbouring 


QL 
powers. 
In practice, to improve the graphical resolution, 
we seldom plot HCN, D]) versus €), nor z({N',D'}) versus 
Yo: Instead the deviations of the estimations from a 


reference line are used in place of H({N,D]) and Ze CINE Do 1) 


that is we use 
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Aw ({IN,D]) = WO CIN,D]) - ae » - b (4.15) 
and 

Az, ({N',D*]) = z€IN',D']) - a'yy te (Phy (4.16) 
where z = ary s tps and. ji = ae » + b are the reference 


lines mentioned. 
Figures 4.12 to 4.15 are graphs of Aus versus €) 
£0r & = 1, 2, 4, 5. ‘The best estamates: of €, are indicated 


in table 4.3. The graphs of Az versus Yo are given in 


2 
figures 4.16 to 4.21 for 2 =) 0, 1). 72, 73, 4,75. he best 


estimates of Yo are indicated in table 4.3. 


It is interesting to note that while the method of 
forming neighbouring powers is more elaborate and precise 


that of taking logarithmic derivatives is much faster. 
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Figure 4.14 Deviations, Au, from a standard line, u=ae,tb, of 
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Figure 4.15 
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Figure 4.16 Deviations, Au, from a standard line, 


z= a! Yo + b', of poles of Padé approxi- 
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Figure 4.17 Deviations, Az, froma standard line, 


Zo eee 4 + b', of poles of Padé approxi- 


OM (u,Z) 
Ou 


1/y3 
versus Yu° 
BR 


20 


_.(9,8) 
(6 2) 


(7, 6) 


(7,5) 
(6,8) 


(6,6) 


(8,6) 


(6,8) 


-10 


Ze ® 0-31740) 


Figure 4.18 Deviations, Az, from a standard line, 


z= a!' ¥5 + b', of poles of Padé approxi- 


mants to 
1/Y5 
Be 2 
{3 oz) versus Y>° 
s) 
bs =u 


\Ze, 2 0: 317401 


“10 


65 


66 


Figure 4.19 Deviations, Az, from a standard line, 


2 aa Y3 + b', of poles of Padé approximants 
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4.4 Transformations 

Already in section 4.2 we have seen that the presence 
of non-physical singularities can prevent the analysis of 
series expansions from proceeding smoothly. First of all 
we wish to find out the locations of these singularities. 
We seek therefore the poles cf central, high degree Padé 
approximants to (€/dz) Log (3°M(u,z)/2n° |, le Figure 4.22 
depicts the case for £ = 4 and figure 4.23, the case for 
& = 5. It appears that there is a pair of singularities 
in the complex z plane making an angle of about 60° with 
the positive real axis. In both these cases these non- 
physical singularities appear nearer to the origin than Zoe 
This unfortunate situation may account for the greater 
error we must associate with our estimates of Y30 Var Yo: 
Moreover, on scrutinizing figure 4.24 we notice the irri- 
tating presence of a negative real singularity in 
(4/dz) Log {97M (w/z) /20 | oy V 

Many ingenious techniques have been developed for 
dealing with non-physical singularities (Gaunt and Guttmann, 
1974). We follow here the method of conformal mapping as 


developed by Betts, Elliott and Ditzian (1971). We rewrite 


the series to be analysed: 


2 17 
0 M(y,Z) Py ) p &*) 2f. (4.17) 
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The most general conformal transformation of 
QR Q ‘ i : : 
[9 M(u,z)/ou ris is achieved by the substitution of 


Cc 


“ = 48 
2 ee ee (4.18) 


where L (L = 17 here) is the maximum degree of the known 
terms in the original series expansion, (4.17). 


The resulting new polynomial is obtained 


Q 


aiM(u,2) ype (4.19) 
hel n= 


For the pair of singularities at an angle of + 60° 


to the real axis in the complex z-plane we notice that the 


transformation: 
- -3 
gos Z/() = 272 ) (4.20) 
has the following effect on points of distance |z| = 1/3 
from the origin. Such points on the rays argument z = 0, 


+ 27/3 are shifted radially nearer to the origin. Points 
on the rays argument z = 7, + 1/3, however, are pushed 
radially farther from the origin. 
As for the case of [a°M(u,z2)/927|_. ], neither 
e 


the technique of logarithmic derivatives nor that of 


neighbouring powers were successful in the estimation of 
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Ee, (see tables 4.1 and 4.3). Most unfortunately, when 

the poles of high degree central Padé approximants to 

(d/du) &n [a°m(u,z)/d271, are examined, a singularity 
appears on the real uw axis ef yu = 0.85. This is rather 
disconcerting and as yet no simple technique has succeeded 
in overcoming this impedance to the analysis of the critical 


behaviour at We = ] for this function. 
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Figure 4.22 Poles of Pade approximants to 
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Figure 4.23 Poles of Padé approximants to 
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4.5 Results of Transformations 
Using the “transformation (4720), ane 0.317401 
is transformed to Ze = Or222720F “The=transtormation (47°20) 
= 4 
was applied to [9 Mey Chae) = and feet Gh eat Hy 
W=U pals 
The resulting estimations of Nien and Ys from residues and 


poles of Padé approximants to 


ea 
ues log * e) Aes! 
dz 2 Ou re 
aie 
and 
d 1 [a®mqu,z) 
— dog 
= -7 5 
dz Zz ou o 
U Ve 


are displayed in figures (4.24) and (4.25) respectively. 

4A 4 5 5 

Next, (a MU, zi7ou | rican iM (ie) ey aie ae 

tae pea 

were analysed by the method of neighbouring powers described 

in section 4.3, after the transformation (4.20) was applied. 
The results are illustrated in figures 4.26 and 4.27 res- 


pectively. We obtain from these figures the improved 
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The transformation (4.20) appears to improve 
appreciably the analyses of the series for [9°M(u,z)/20) 
and [97M (u,2)/307] . There is a marked improvement in : 
precision as shown by the better confidence limits we attach 


to ag and Yes The best estimates of Ep and Yo in the overall 


analyses are given in table 4.4. 
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4.6 Calculation o£ Critical Amplitudes 
We have so far analysed the asymptotic behaviour 
of (9°M/au") and (a°M/az") . The expected behaviour 
ie c 


of the most singular parts of these functions may be 


expressed as: 


CC VENC ey Tome faley, nm (4.21) 
UFU sg c 
ve g = z smote 
(0 M/dz CE % Ey (al Lu) Q Cae) 


Theories of critical phenomena are usually 
expressed in terms of the reduced temperature t = 1 - T/T oe 
and the reduced magnetic field h = MH/KT for a maqnetic 
system. The expected critical behaviour in terms of t and h 


should then read: 
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: 1 ory R a s Le 
ot Ch (=2) (Coys log z) Q (4A°25) 


We note that log ZO is negative, so (-log Zo) is positive. 


Upon equating (4.22) to (4.24) we get: 


eo) Am omens ly A R 
(Ey) Ne (E,) had oy le (Z Log Z 


res 
II 


as € yy 
(Ey / Jee) (ZO log ZO) (4.26) 


The resulting values of Eo and Cy estimated are given in 


table 4.4. 


et iw: 
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We have seen that the series (2 = 0, 1, 2, 3, 4, 


17 

L 
[a'M(y,z)/au-] = J 
rN n=0 


has its leading non-zero coefficient, Be aie Eso 


for £ > 0. For £ = 1, our truncated series in (4.27) has 

15 non-zero terms. Finally, when £ = 5, the series, for 

aman] has decreased to 11 non-zero terms. We know 
Cc 


that the important terms are the higher powers of z, for 
example terms in z! where p > 17, Also in general, the 
longer the series the more accurate the estimations made 
from the series. To obtain terms involving z?, where 
p > 17, we would need to have higher L, (2) polynomials. 
Since we only know exactly the L (2) polynomials up to 
s = 23, we cannot, for the moment, extend the series in 
(4.27) beyond el However, is there a way of increasing 
the string of non-zero terms in (4.27)? 

One method of retaining as many non-zero terms in 
Z as possible is to study instead the series expansions 
representing the asymptotic behaviour of (u9/34) M| . 


Recall equation (3.11): 
M(u,z) =1-2 } s L,(z) ue 


uw OM/ou = =— 2 ) sb. (z)ue 


(2) n 
bo zZ (4.27) 
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In general 


ole 23 
Cima MU e) ie = 21) 


L (Zz) uo (4.28) 
s=1 


In particular, on the coexistence curve: 


23 
One. L+1 
(u x)" M(y,z)|,_. =-2 } 8 L. (z) 
ou HE s=1 
(4.29) 
Ly 
= sae) ae 
sai 2 
ry) 4(2) op 
The table of coefficients for - 2 } eo Dowd aoe Oe Lea es 
=0 


3, 4, 5, are given in appendix 6. 


For 2 = 0, 1, the series representing (H9/3u) MIL, 


and the series representing a *Msaur | are the same. 


y= 
rs! 
For £ = 2, the series representing (18/8u) 7m] has 15 non- 
c 
zero terms while the series representing 37 M/ aye | = has 14 
ce 


non-zero terms. The truncated series 
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(2) my be 
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was analysed by computing poles and residues of Padé approxi- 


mants to: 


sas SII lin zZ : (4.31) 
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A 3 by 4 Padé (approximant) table was set up and the 
computed poles and residues are listed in appendix 7. 


Y>5 estimated by this method was found to be 
Y¥oe= 22988 2° 0, 10 (4.32) 


Next, poles and residues of Padé approximants to 


Dey S NS 
| ) ee : 4 were computed. By this method, for 
n=3 


a range of values of Y3 such thar. 27.90 "< Y5 Se OO rae 
was found that 


Yio ee 0704 (4,33) 
Z 


In the case of 2 = 3, an analysis of the poles and 


residues to 


: Jky) 
— log ) eae ae ; (4. 34) 
n=3 
yielded: 
¥3 ae eee O08) (42.35) 


It was suspected from the large uncertainty in (4.35) that 
troublesome non-physical singularities were hampering the 
analysis of the series. A natural step to take is to 


examine the locations of the singularities in the complex 
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z-plane of central high degree Padé approximants to 

(4034) ,2 fOr 82 =" 2, 3; 4, 5. However for << ="3,. 47-5, 

it was found that the series in (4.29) cannot be readily 
analysed by the usual techniques of series analysis whereas 
the series in (4.27) can. We have therefore chosen to 


concentrate on the analysis of the series in (4.27). 
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4.8 Dropping Off Terms 


We have also attempted to investigate the effect 
of dropping off leading non-zero terms in the field deriva- 
tives of the magnetization on the coexistence curve. For 
example the known non-zero coefficients in the truncated 
series expansion for 

17 


(97M /2u7) = y see 3m : 
c n=0 


are Bi*? to bees From the voles and residues of Padé 


approximants to the logarithmic derivatives of the modified 


series: 
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We found that the estimated Y5 fluctuated between 2.95 
and 2.55. We realise that the important terms are those 
involving the higher powers of z. If this is the case then 
dropping off the terms involving the smaller powers of z 
should not affect adversely the asymptotic behaviour of 
the series at all. 

This procedure of analysis was repeated with the 
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Beginning with the spin 1/2 Ising model on the 
hydrogen peroxide lattice, the low temperature (high field) 
series expansion for the Gibbs free energy has been 
obtained from the configurational data of Betts, Elliott 
and Sykes (1974). The magnetization and its first five 
field derivatives on the coexistence curve, and the magneti- 
zation and its first five derivatives with respect to 
Z = exp(-2J/kT) on the critical isotherm have been derived. 
Next, using the techniques of taking logarithmic derivatives 
and neighbouring powers, the series for (B°M/9u") and 
eo have been analysed. Where the ee has 


c 
been impeded by the occurrence of non-physical singularities 


(a’M/az~) 


near the physical singularity, transformations have been 


sought to facilitate and enhance the analysis. From the 
4 
Psael 
fear 


residues of Padé approximants to (a*M/au and 


L/e 


Q Q ve 
(d M/oz ae 


, the critical amplitudes Cy and E, (Cina a eu 
and 4.28) have been computed. Finally, the best estimates 
of Yo and Ey (2 = 0, 1, 2, 3, 4, 5) have been tabulated in 
table 4.5 together with those predicted by scaling theory 
for the purpose of comparison. 


If the series for (Q°M/au") a, are examined, we find 
‘e 


that for approximating M(z) (i.e. the case & = 0) we have a 


series of 18 terms (i.e. ) b Zz At the other extreme, 


for 2 = 5, we are approximating (9°M/9u") by a series of 
7 (2 

only eleven terms, (i.e. )} Be z), This means that when 

n=7 

we form Padé approximants, to the logarithmic derivative of 
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(39°M/ay°) cs or to Bey Re als , we are restricted to 
UH UU, 

a smaller Padé vapte than in the case of 2 = 0. This may 

explain the general trend towards a smaller degree of 

precision (i.e. wider confidence limits) as 2 increases. 

Nevertheless, as has been pointed out by Betts and Chan 

(1974), all the critical exponents ic Citra Oe i aa. Son okt) 

are in good agreement with static scaling predictions where 

it has been assumed that o = 1/8 and y = 5/4. 

From table 4.4, the Yo values agree within confidence 
limits with scaling predictions for 2 = 0, 2, 3, 4. In the 
cases of 2 = 1, and 5, the disagreement with scaling predic- 
tions is never greater than 4%. The confidence limits tends 
to increase with increasing 2 which is a reason why investi- 
gations were confined to 2 < 5. On the whole the critical 
exponents, Yor along the coexistence curve are in good 
agreement with scaling predictions. 

In the case of €, on the critical isotherm, Betts 


£ 
and Filipow (1972) have investigated €, for a number of two 


gL 
dimensional lattices: the honeycomb, the square and the 
triangular lattices, using the configurational data of Sykes, 
Gaunt, Mattingly, Essam and Elliott (1973). They found it 
impossible to estimate Ey at all. They obtained quite pre- 
cisely, estimates of Egr &5 and Eger which agree closely with 
scaling predictions. Their values for Ey and €, were not so 


precise but nevertheless the agreement with scaling predic- 


tions is still satisfactory. For the three dimensional 
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hydrogen peroxide lattice it was not possible to estimate 


E3 at all. Also, while we could estimate €, by the method 


0 
of logarithmic derivatives the method of neighbouring 


powers failed completely to yield an estimate of ¢€ For 


0° 
&£ = 0, 4, 5 in the case of the hydrogen peroxide lattice, 
the agreement with the corresponding e,'s of scaling is 
excellent. For a and E, the agreement is fair. It is 
clear that there is some qualitative similarity between the 
two dimensional cases and the three dimensional hydrogen 
peroxide lattice case. Most of the critical exponents agree 
very well with scaling. Some agree only satisfactorily, 
while a few are not even estimatable. 

On the whole we conclude that the results obtained 
in the analyses of the critical exponents, Yo and Eo for 
the hydrogen peroxide lattice support static scaling for 
the three dimensional Ising model. To obtain even more 
precise results we have to calculate more terms in the low 
temperature series expansions, especially for the field 
derivatives of the magnetization on the coexistence curve. 
A subsequent extension of this work is to test lattice- 
lattice scaling, or universality theory, using the critical 


amplitude data obtained here for the hydrogen peroxide lattice 


and elsewhere for other lattices. 
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APPENDIX 6 


Poles and residues of Padé approximants to 


= log : eel ali 
n=3 
POLE Residue 

PADE (6,7) 0. 339457840990814 —0'56 7965266465912) 
PADE (7,6) 0232 3570550845277 -0.334813316940444D 
PADE (5,7) 0.295750439416026 -0.185440516806951D 
PADE (6,6) 0. 316746646694732 -0.294907918467118D 
PADE (7,5) 0. 2'36112700532900 i=0.932988832077709D 
PADE (4,7) Of G0S28092Z5 282578 (25900271 LO458L6) 
PADE (5,6) Oes0L II TeS 4662018 =0.2033382097354968D 
PADE (6,5) Us32/2031526 74791 -0.388979824836032D 
PADE (4,6) 0.298259784449 363 =021920920342393890D 
PADE (5,5) 0.289504969941134 -O0159349722735093D 
PADE (6,4) 0. 34961267409 3252 =-0.6/3113035637042D 
PADE (3,6) On 2917653713437 21 =0.170975483072717D 
PADE (4,5) 0.249918303485047 - 0593473120997 49 70D 
PADE (5,4) 0.276290868518204 -0.119025228233482D 
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